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Abstract: We consider a holographic QCD model for light mesons beyond the leading 
order in the context of 5-dim gauged linear sigma model on the interval in the AdSs space. 
We include two dimension-6 operators in addition to the canonical bulk kinetic terms, and 
study chiral dynamics of vr, p, ai and some of their KK modes. As novel features of dim-6 
operators, we get non- vanishing Br(ai — > tt'j), the electromagnetic form factor and the 
charge radius of a charged pion, which improve the leading order results significantly and 
agree well with the experimental results. 
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1. Introduction 

To understand the dynamics of low lying hadrons from underlying QCD has been a long 
standing problem in theoretical physics. In chiral Lagrangian approaches, it has been 
known for some time that the low energy dynamics of pions, vector mesons p and axial 
vector mesons oi are well described by the gauged linear sigma model (or its nonlinear 
version) with massive Yang-Mills gauge filds. The model Lagrangian up to dimension-6 
operators is given by ^ 



-^MassiveYM — Tr 



1 



1 



+ Tr 



+ AiTr($t$)2 + ^2 T¥($t$) + (A3det($) + H.c 



(1.1) 



""^We ignore the Wess-Zumino-Witten term in this work. 
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The role of the higher dimensional operators in light hadron dynamics, especially k and 
^ terms, were studied in the framework of the gauged linear sigma model in 4D |2|. 
Although the above Lagrangian is quite successful in describing the n — p — ai system, 
it has a conceptual drawback in that we need to give gauge boson masses ttt-q by hand. 
If we put rriQ = 0, global chiral symmetry becomes local symmetry, which is not a true 
symmetry of real QCD, and we end up with massless p and ai, which is phenomenologically 
disastrous. Therefore we have to put uiq ^ and have to impose chiral symmetry only as a 
global symmetry. However, if chiral symmetry is only a global symmetry, then there is no 
compelling reason to introduce gauge covariant derivative, hence no reason for the minimal 
coupling between hadrons and (axial) vector mesons, and thus universality of the P—V — V 
couplings. Since the universality seems to hold to a good approximation, it is tempting to 
implement global chiral symmetry to local symmetry. This have been remained a problem 
in chiral dynamics approach to the low lying hadrons. 

Recently, there have been many interesting and successful attempts to understand 
hadron physics in the context of AdS/CFT correspondence The properties of hardrons 
and the hadron physics phenomenology are studied in various approaches |5|, |6|, |^, |8|, |9|, 
1^, |ll], [l^, 0] which are inspired by the AdS/CFT correspondence. 

One may start with some stringy setup that may reproduce certain aspects of non- 
perturbative QCD. The most successful approach so far seems, arguably, the works by 
Sakai and Sugkimoto IH, and follow-up papers The model by Sakai and Sugimoto 

has nice features, but also some drawbacks. They show that chiral dynamics of vr, p and 
ai system can be well reproduced by studying the Nj D8-branes in the background of 
D4-branes. Also the Wess-Zumino-Witten (WZW) term is derived from the 5-dim Chern- 
Simon (CS) term. On the other hand, there are spurious SO{b) symmetry from 5^, which 
is not a true symmetry of real QCD. And it is not easy to accommodate nonzero quark 
masses, namely nonzero pion mass. Finally the pion and its radial excitation comes from 
different 5-dim fields, which is not easy to understand within the quark model. Despite 
numerous remarkable successes of Sakai-Sugimito model, there is an ample room for further 
improvement. 

Independent of the stringy approach, a gravity dual model of the gauged linear sigma 
model was proposed to describe the chiral dynamics of light hadrons [P, 0]. This approach 
is often called the bottom-up approach, where one starts from QCD and then tries to 
construct its five-dimensional holographic dual model, AdS/QCD. Following the AdS/CFT 
correspondence, it is assumed that there are bulk fields that couple to the 4-dimensional 
QCD operators. For example, there are bulk gauge fields Lm and Rm that couple to the 
QCD operators ji = QlI^Ql and Jr = qn^^QRi which are flavor currents. 

Quality of the overall fit to the meson properties in the models of Ref. |^ is at the 
level of ~ 30%, which is quite remarkable, considering the simplicity of the model. However 
it predicts B{ai ttj) = and too small charge radius of a charged pion. The Lagrangian 
in Ref. is the leading order one, since it contains only the bulk kinetic terms for the 
bulk gauge fields and scalar fields. In order to improve the predictions for B{ai — > ir^) and 
the charge radius of a charged pion, we have to go beyond the leading order Lagrangian. 

In this paper, we construct an AdSs dual model of the gauged linear sigma model 
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with dimension-6 operators, motivated by the recently developed AdS/QCD model ||6|, 0. 
To this end we incorporate higher dimensional operators, especially two dim-6 terms, into 
the AdS/QCD model. ^ In this work, we only consider the vector, axial- vector and 
pseudoscalar sectors, as a first step of our study. Interestingly enough, we find that the 
aforementioned problem of giving gauge boson masses tuq is no longer present, since one 
can give masses of the vector and axial vector mesons, by projecting out the zero modes 
by choosing suitable boundary conditions. Degeneracy between the vector and the axial 
vector mesons will be lifted by the conventional Higgs mechanism. Still there remain 
physical pions. 

Naively, these new operators will have nontrivial effects on the interaction vertex, as 
well as mass spectra and decay constants. We expect that they may contribute to the 
Br(ai 717), which is zero in the orig inal AdS/QCD model |, 0. We also study the 
phenomenology of p ^ vrvr and ai — > pvr, the branching ratios and D/S wave amplitude 
ratio in the latter channel. By introducing photon as an external field, we study the pion 
electromagnetic form factor, and calculate the pion electromagnetic charge radius, which 
agrees with the experimental results in our numerical study. 

This paper is organized as follows. In section 2, we define the Lagrangian of our model 
with two dim-6 operators in AdSs. In section 3, we study the mass spectra and decay 
constants in vector, axial-vector and pseudoscalar sectors. We also present the interaction 
vertex and phenomenology of ai — > pir, p — > vrvr, oi — > irj channels, and calculate the pion 
charge radius. We derive the relevant chiral coefficients in section 4, and give our numerical 
results in section 5. The conclusions are drawn in section 6. 

2. Gauged Linear Sigma Model in the AdSs space 

The Lagrangian of the holographic QCD model |^, ^ defined in a slice of AdSs is given by 



/:dim-4 ^ Tr 



1 T tMN 1 7? TjMN 



(2.1) 



where M| = -S/L^ from AdS/CFT correspondence §], Dm^ = + iLm^ - i^Ru, 
Lm = L%^t''/2 with being the Pauh matrix, and M, iV = 0, 1, 2, 3, 5(or z ). We define 
$ = 5e^^/K^) with {S) = v{z). Under SU(2)^, S and P transfor m as singlet and triplet, 
respectively. The AdSs space is characterized in the conformally flat metric with a warp 
factor a{z) = L/ z: 

ds^ = a'^{z){dx''dx^ - dz^). (2.2) 

The scale L is the curvature of the 5-dimensional AdS space. In this model, the AdSs space 
is compactified such that Lq < z < Li, where Lq ^ is an ultra-violet (UV) cutoff and 
Li is an infrared (IR) cutoff. Solving the equation of motion for S, we obtain ||^] 

{S) = v{z) = ciz + C2Z^ (2.3) 



^For a brief report on the present work, see Ref |16l 
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with the integration constants ci^2, 



ci 



MgLl 



LL,{Ll-Liy LL,{Ll 
Here we adopted the fohowing boundary conditions 



i - MM 



(2.4) 



L 



i = Lv 



Lo 



(2.5) 



where Mg is the current quark mass matrix, which breaks chiral symmetry exphcitly, and 
is related to {qq), which breaks chiral symmetry spontaneously. The value of Li is fixed 
by the rho-meson mass: 1/Li ~ 320 MeV There may be several ways to improve the 

model given above, though several observables obtained from the model are in agreement 
with experiments. One immediate extension of the model is to see corrections from various 
sources: trilinear or quartic interactions among the vector fields, 5D loop corrections, 
higher dimensional operators and back-reactions on the metric due to condensates |12]. In 
the present work, we consider corrections to the model from higher dimensional operators, 
though to be consistent we have to treat all those corrections at the same time. We 
note that a part of large Nc corrections through meson-loop contributions are discussed in 



Ref. lU. 

Now we introduce higher dimensional operators in the model Lagrangian in Eq. (p.lj). 
In principle, we can include infinite tower of higher dimensional operators, but for simplicity 
we consider only dimension-6 operators in the chiral limit. Note here that we have the 
following mass dimensions for a scalar field $ and vector fields Lm and Rm- 



dim(<I)) = dim(LAf) = dim(i?M) = 1 
The Lagrangian with dimension-6 operators reads 



(2.6) 



dim— 6 



V9M5 Tr 



C 



Mi 



where k and C are constants that will be fixed later. 
There are more dimension-6 operators, such as 



(2.7) 



^5 



V9M5 Tr 



M 



(2.8) 



However these terms are 0{p^) after chiral symmetry breaking, whereas the k and C, terms 
are 0{p'^) after chiral symmetry breaking. Therefore we keep only those dimension-6 
terms that reduce to O(p^) after chiral symmetry breaking. We note that the corrections 
to physical observables from the second dim-6 operator Tt[Lmn^R^^^ ^'^] in Eq. (2.7) and 
the operator in Eq. (2.8) have been discussed in Ref. |jl9|. 
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3. Vector, axial-vector and pseudoscalar sectors 



3.1 Relevant parts of the Lagrangian 

In this section, we work in the chiral hmit. Then v{z) is proportional to 1, v{z) ~ Cf^l- 
The vector and axial gauge bosons are defined by 

1 



1 

Am — —^{Lm - Rm) 



(3.1) 



In order to cancel the mixing terms of V^, (/Li as 4D Lorentz index, 0,1,2,3) and Vz 
Az, P, we add gauge fixing terms 



GF 

A 
GF 



M50 



Tr 
Tr 



a 



2C 



d^iav^Vz 



df.A" - ^ (d^iaAz) + V2a3t;P 



In the unitary gauge, (,v,A — > 00, we have the following relation between Az and P, 



(3.2) 



2a/2k 4k 2C 

V2a^vP + 55(0^5) + ^79-^5 (a(55i')-P) + — 2aw(55t')A^ + -^55(au2A^) = 0, (3.3) 



Mi 



which is identical to the leading order relation when k = C = 0. 

The quadratic terms for vector, axial-vector and pseudoscalar are given by, after inte- 
gration by parts. 



Ms 
2 



Ca = ^aTrl^A^(^d'Za - a-'d,aZ,d, + 2a\^ - ^v{d,v)d,^ A^^ , 



^aTr|( 



,-^-^(-2aV)^^. + ..^ 



P 



+a{d,Az)^ + a\d^Pf + ^:^a{d,v){d^Az){d>^P) + ^^av\d^Azf'^ . (3.4) 



with Zy = 1 — and Za = I + ^jr- The boundary terms are 



^boundary = MgoTr (^^^Z^^s^^M + ^''^a^S^p 



-A„,^^^A,, - '^^{d,v)A^d''P - ^v'^A^d'^A, 



z=L\ 



. (3.5) 



z=Lq 
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We choose the following boundary conditions to cancel the IR-boundary terms, 



z=L\ 



0, 



z=L\ 



A, 



^z + ^i^{d^v)P+^v^A, 



= 0, 



z=L\ 



(3.6) 
(3.7) 



z=L\ 



and the UV-boundary condition will be specified later. 

We also calculated VAP, VPP and four-pion interaction vertices 

V2i 



2iK ^ 
-Tl-aTr 
M5 



jO-vap = ^MsoIV A''[d5V^,A,] - {d5A^')[V^, A,] + V2a\A^'[V^,P] 

id,v)id5A^')[V^,P] - V2vid5v)A''[V^,A,] + V2v^A''[d5V^,A, 
+vA''[d5V^„ d^P] + vAi'lV^,, d^P] 

Af'^[V^,,P] - 2{d5A'')[d,Vf„P] + V2vA^[d5V^,A,] 



M5 



+V2v{d5A'')[V^,A, 



(3.8) 



■\/2iK 



Ms 



V2iC 

Ms 
a^Ms 



V^[A,,di,A,]+a^V^[P,d^P] 



lv'"'[d^P,d,P] + {d5V'')[d^P,d5P] - V2v{d5V'')[A„d^P] 



+V2{d5v)V>'[A„d^P] - V2{d5v)V''[d^A„P] 



avTr 
TV 



- vVf'iA,, d^A,] + V2id5Vf')[P, d^A, 



(3.9) 



{d.pyp' - {d^p)p 



(di'A,)(d^P)A,P-\{d„A,)P 

V2 a{d5v)K 



3 u^Ms 
\f2aK 



-Tr 



(a'^A,)(5^p)pp - {d^A,)p{d^p)p 



vMr, 
< nn 

-T^Tr 

Ms 



-Tr 



{d'^A,){d^p)p{d5P) - {d^^A,){d^p){d^p)p 



{d^A.fp'' - ( {d^,A,)P 



2n 



(3.10) 



3.2 Two-point correlation functions 

We calculate the two-point correlation functions for vector and axial-vector with respect 
to the UV boundary external source fields 1;^ and a^, which couple to the vector and 
axial-vector currents operators, respectively, 

V^\z=Lo = V^, A^\z=Lo = «M- (3.11) 
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From the AdS/CFT correspondence, in order to calculate the current-current correlation 
function in the strongly coupled CFT side, we can do it in the weakly interacting AdS side 
instead. Then the effective Lagrangian in momentum space in term of the correlators is 



(3.12) 



with Uy'^^p^) = {g'^'^ — p^^p^ /p^) Hy^Aip^)- We solve the equations of motion for the vector 
and axial- vector field derived from Eq.( |3.4j) with the boundary conditions (|3.6| ) and ( 3.11] ) 
and calculate the two-point correlation function 



n(p^) 



(3.13) 



z=Lo^O 



where f{z) is the solution of differential equation. With dim-6 operators, we cannot cal- 
culate the two-point correlation function n(p^) analytically, instead, we do it numerically. 

For asymptotically large momentum p^L^ ^ 1, we can expand the 2-point functions 
in powers of 1/p^, and get 



ny,A(p') = 



M5L 



Inp^Lp -I- c, 



with 



4 



192C 



4 



/I6M5L 192C 



1 



384k 



(3.14) 



(3.15) 



which agree with the results in Ref. for k = = 0. It is worthwhile to calculate the 
left-right correlator II^j^ = Ily — U^, in the large momentum limit, we have 



£6 

A 



P 



(3.16) 



with cq = cv — CA, where the experiment value cg = — 47ras(gg)^ ~ —1.3 x lO^'^GeV^ 
obtained from Ref. [^] . We remark here that the vector correlator obtained in the present 
work and in the hard wall model |^] has no 1/p^ compared to the results from operator 
product expansion (OPE) |21]. In the chiral limit, the coefficient of 1/p^ term is due to 
the gluon condensate ||2l|. In the hard wall model adopted in the present work, however, 
the metric is just a pure AdS with no gluon condensate included, and the model has no 5D 
bulk scalar field that couples to iv{G ^j^yG^^) at the boundary, where G^y is the gluon field 
strength tensor. Therefore, the vector and axial-vector correlators in the hard wall model 
do not contain 1/p^ term, as it should be. To have 1/p^ in Eq. (3.14), we have to consider 



a deformed AdS background [22| due to the back-reaction of the gluon condensate. 

In the large Nc limit, the above correlators can be written as the sum in terms of the 
resonance masses and decay constants. 



nA(/ 



nv(/) 



p 



n ■ 



/a. 



Ml 



+ fl 



f2 



Ml ■ 



(3.17) 



(3.18) 
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Then the vector and axial meson masses are determined as the poles of their corresponding 
correlators, and the decay constants are related with the residue, 



lim (p^ 

nA(o). 



m'p^,^)Ilv,A{p)/P 



(3.19) 
(3.20) 



4. Interactions and phenomenology 
4.1 KK decompositions 

Now we study hadronic observables such as decay widths and form factors using our model 
given in Eq. ( |2.1| ) and Eg. ( |2.7| ). Primarily we investigate how those dimension-6 operators 
in Eq. (2.7) affect the results obtained with only interactions in Eq. ( p.!] ). To this end, we 



first Kaluza-Klein (KK) decompose the vector field as V^(x, z) = ^^^^^ 
and also for the axial- vector and pseudoscalar fields, where we omit the superscript index 
(n) when we consider the lowest KK mode. The first resonances of the vector, axial-vector 
and pseudoscalar fields are associated with p, ai and vr respectively. The equations of mo- 
tion for the vector, axial- vector and pseudoscalar fields are easily read off from Eq. (^.-j|). 
To cancel the boundary terms, in addition to the IR boundary conditions given in Eq. ( |3.(j ) 
and Eq. (|3.7D , we impose the following UV boundary conditions 



V„ 



0, A, 



z=Lo 



0, 



z=Lo 



0. 



(4.1) 



z=Lo 



We obtain the wave function and mass spectra of various fields numerically with the nor- 
malization conditions: 



dz^Zviz)fy^\z)f^\z) =5mn, 
dz^Za{z)f^^\z)f^^\z)=Smn, 





dz— 


r 




dz— 




r 


dz— 



4V2k 



1. (4.2) 



4.2 p TTvr 

The pmr vertex can be expressed as 
with the couplings 

+^ ( - ^^^fv)fp{d^fp) - V2v{d^fv)fAjp + 2V2{d5v)fvfAjp^ 



9pn 



fvfl+a^fvfp 
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( - vfvfl + V2{d,fv)fAjp^ 



(4.4) 



fpnn = I dz ^l^ afvfp (4.5) 

We also calculate the decay width r(/> — > vrvr), which includes the non-minimal coupling 
fp-Km even though its contribution is numerically small. 

4.3 Electromagnetic form factor of a charged pion 

Before we study the electromagnetic form factor of a charged pion, we introduce the photon 
as an external gauge field and rewrite the bulk vector field decomposition as 

^ oo 

V,{x,z) = eF,{x)T, + ^=Y^V^-){x)f^^\z) , (4.6) 

^ n=l 

with Ts = o"3/\/2, where a is the Pauli matrix, and e is identified with the physical electron 
charge at chiral symmetry breaking scale. To treat photon and p on the same footing, 
we introduce friz) = 1 as the fifth dimension profile for photon. The advantage of our 
treatment of photon as external field, compared with the treatment of photon as the electro- 
magnetic subgroup of SU {3)v 0, is that we don't need to worry about the KK excitations 
of the photon, as well as the mixing between photon KK excitations and KK excitations. 

We consider the electromagnetic form factors of pions. In additional to the usual 
structure of contact ^tttt interaction Tr(F^[A5, df^A^]), we also have non-minimal structure 
Tr (F'^'^ldfj^A^, di^A^]), which comes from the dim-6 k term. And we also find g^n-K = e after 
comparing with the pion normalization condition, eq.( [4.2| ). From the kinetic term ( |3.4| ) 
and the vector KK decomposition ( [4. 61) , we can derive the kinetic mixing of 7 and p, 



2 

with 

r-L 

The electromagnetic form factor of pion can be calculated as 



M5 /■'"I 

3^p^~fT^ / dzaZyfviz). (4.8) 



F{q') = 1 - - -P^gp... (4.9) 



In small momentum limit, it can also be expressed as 



F{q^) = l + Klq' + 0{q^), (4.10) 



with the pion charge radius calculated as 



r! = 6 



(4.11) 



Our vector meson dominance (VMD) is difi'erent from the usual VMDs as discussed in 
ref. 1 24], where we have an additional non-minimal 'jTrn contact interaction, Ti {F^'^[dfj^A^, duA^,]). 
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4.4 fli —>■ pTT 

We first consider the process oi — > pvr. Applying the KK-decomposition to Cvap in Eq. 
(13.81), we obtain 



(4.12) 



+^g3a^pnTT{A^'''[V^,,A,]) 
with the coefficients gia^pn (i = 1,2,3) 



9lai 



pn 



dz 



2k 



-^[fAidMfA. - id5fA)fvfA. + V2a\fAfvfp 



-^\{d^v){d-,fA)fvfp - V2v{d5v)fAfvfA. 
+V2av^fA{d5fv)fA, + avfA{d5fv){d5fp. 

'-'^ '{dMfvfA. + vfAid5fv)fA. - V2{d5fA){d5fv)fp 



92ai 



pn 



az- 



gSaipn 



az- 



avfAfvfp 



av^{d^fA)fvfA, 



(4.13) 
(4.14) 
(4.15) 



With the interaction vertex above, it is straightforward to derive the amphtude of the 
process, which can be written as 

The S/D wave amphtudes are defined as in Ref. [0] 
{p{ks,)7r{-k)\H\ai{0sa,)) = loo(f^fc) + J] C(211; mz.5pS,jy2m Jf^fe) 



faipTrgfiu + gaipTrPiT nPnu 



with 



J aipn 



3m, 



{Ep + 2mp)fa^pT, - k'^rua^ga^ 



D 



3m, 



{Ep - mp)fa^p^ - k'^nia^g 



'a\p-K 



And also the decay width of oi — > /ovr is 



r(ai pn) 



Pc 



ai 



1 / E 



m 



Q ./aipTT 2 1 ^ JaipTvn- ^ Pcgaipi^ 



2-1 



(4.16) 



(4.17) 



4.5 ai — > TT'j 

In this subsection, we study the process ai — > TT'j, With the help of vector KK decom- 
position Eq. (|4.6| ), we have similar results as oi — > pir. We have verified that the gauge 
non-invariant term of structure Tr { A^^[Fp, A z]) is cancelled out, when we impose the rela- 
tion between Az and P, e.g., Eq. ( |3.3[ ), the boundary condition Eq. ( ^ ) and d^fp = 0. 
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4.6 Numerical results 

In this subsection, we present the numerical results of various hadronic obsevables and 
chiral coefficients discussed previously. We use to fit the four parameters Li, ^, k, ( 
from nip, rua^, D/S ratio, T{p — > vrvr) in case B and nip, niai, ^(P ~^ ^^)> /vr case C. 
Our results are summarized in Table |l|, 0, and Figure ||. As a comparison, we also give 
Da Rold and Pomarol's results in case A. 

In both cases B and C, T(ai tt'j) is non- vanishing, but small (less than lOOKeV), 
while r(ai pn) is a little small in case B, but consistent with experimental measurement 
in case C. We have checked that the dominant contribution to T(ai pn) comes from 
the leading order structure Tr(j4'^[y^, ylz]). However, Tic{A>^[Fp, A^]) term is not gauge 
invariant and cancelled out for oi ttj channel. Then only dim-6 k and C terms contribute 
to the above process. This is different from usual 4D models with large r(ai t^j), where 
the ratio between T{ai — > ir^) and T{ai pn) is roughly /gj^^^, and only a single type 
of operator Tr(^'^'^[V^i,, vr]) contributes to both channels 



case 


Li 


K (10-6) 


■nip 


ITlai 


r(p — !■ nn) 


r(ai nj) 


r(ai pn) 






C (10-^) 


fp 


fax 




f7r(fm) 


D/S ratio 


expected 






775.8 ±0.5 


1230 ±40 


146.4 ± 1.5 


0.640 ±0.246 


250 ~ 600 


86.4 ±9.7 












0.672 ±0.008 


-0.108 ±0.016 


A 


3.125 


0. 


[769.6] 


[1253] 


95.4 


0. 


295.5 


85.0 


4.0 


0. 


138 


163 


4.8 


0.585 


-0.055 


B 


2.836 


-5.930 


[775.8] 


[1230] 


[146.5] 


0.088 


165.3 


71.9 


2.56 


-39.72 


144 


182 


5.8 


0.654 


[-0.094] 


C 


3.102 


-16.03 


[775.8] 


[1246] 


[146.4] 


0.042 


409.8 


[78.7] 


4.010 


0.09188 


140 


172 


5.6 


0.640 


-0.027 



Table 1: Various hadronic observables obtained in the present work. The unit of masses, decay 
constants and decay widths is MeV. The inputs for each case are shown in the brackets. 



The pion charge radius r^r agrees with the experiment in both case B and C. Pion 
decay constant /,r is a bit small in case B, while the D/S ratio of ai pn is small in case 
C, compared with experiment. As in other 5D models, the KSRF relation gj^^^/m^p = c/ f"^ 



with c = 1/2 [23] is not satisfied very well. In both case B and C, c is roughly 0.3, which 
means the complete vector meson dominance of order 0{p'^) four-pion interaction, with 
the higher p resonace and scalar exchange, and contact four-pion interaction contribution 
below ~ 10%. 





case A 


case B 


case C 


expected 


^6 


0. 


0.0008 


0.0000 


-0.0005 


„A 
^6 


0.0014 


0.0000 


0.0006 


0.0008 


C6 


-0.0014 


0.0008 


-0.0006 


-0.0013 



Table 2: OPE coefficients Cg', c^, and cq in unit GeV^. 



The pion form factor F{(f') as a function of is plotted in Figure ||. We find the form 
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Figure 1: Pion form factor F{Q'^) as a function of q^. The white circles are data from CERN p5| , 
square from DESY , triangle from DESY |^ . black circle from Jlab |^ , and black square from 
Jlab |§. 



factor has better behavior in case B and C for large value of momentum than that in case 
A. 

The OPE coefficients , Cq, and cq are presented in Table |. The individual coeffi- 
cients and do not agree very well with the expected value in all three cases, while 
the coefficient of left-right correlator agrees with the expected value in case A. However, 
we note that the OPE is also sensitive to the deformation of the AdS metric Q . Consider- 
ing the OPE behavior, it may be worth remarking that, although the spirit of bottom-up 
AdS/QCD models has been to match the theory in the UV and then compare with the 
physical observables in the IR, it is not surprising that the best fit to data would arise 
from a model that disagrees with the precise UV behavior of QCD, where the model is not 
expected to be valid. 



5. Chiral Lagrangian for pseudoscalars up to 0(|)^) 

Before we discuss the 0{p'^) chiral Lagrangian, we consider the vector field p effective 
Lagrangian [|^, 

^ e5^pT¥[V/f ] + ^/p,,/2TV[Vn'^nn (5.1) 



with p transforming as gauge field of SU(2)v and the notation of P^, Z^'^, the same as 
in ref. [30|. The coefficients in the effective Lagrangian are determined by matching with 
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our theory with dim-6 operators. The ©(p^) chiral Lagrangian for the pions is given in 
ref m, 

Ca = LiT?[D^U^ D^'U] + L2Tt[D^U^ DuU]Tt[D^'U^ D^U] + L3TT[Df,U^ D^'UD^U^ D'^U] 
+ UTr[D^U^D''U]Tr[U^x + X^U] + L^TrlD^U^Dt'UiU^x + X^U)] 
+ L6Tr2[[/tx + x^U] + LrTr^ph - x^U] + LsTr[x^Ux^U + U^U^x] 
- iLgTv[F^''D^UD,U^ + F^'D^U^D.U] + LioTt[U^ F^'UFl^,] . (5.2) 

In the present, we do not discuss scalar and pseudoscalar resonances contribution to 
-^3,4,5,6,7,81 ^-nd only study the vector and axial resonances contribution to -/^i, 2, 3,9,10- After 
Integrating out the vector rho meson, we obtain the following chiral coefficients, 

T J TT 2 J IT 



OpTVTT ~ -7-^9 piTTT f p-KTT , ^2 — 2Ll , L3 — — 6L1 



T J TV 2 I J TV £ J TV /r q\ 

Lio can be calculated from the two-point correlators of vector and axial, Ily^^, 

Lio = ^[n'^(o)-nV(o)], (5.4) 

where the derivative is over p^. 

We also calculate the electromagnetic mass difference of the pions from the operator 



3ap 



dp^(nA-ny). (5.5) 



case 


Li 


L2 


is 






TO^+ — m^o (McV) 


exp 


0.4 ±0.3 


1.4 ±0.3 


-3.5 ± 1.1 


6.9 ±0.7 


-5.5 ±0.7 


4.6 


A 


0.43 


0.86 


-2.6 


5.1 


-5.5 


3.4 


B 


0.32 


0.65 


-1.9 


4.0 


-5.0 


1.5 


C 


0.46 


0.93 


-2.8 


5.3 


-5.1 


2.9 



Table 3: The chiral coefficients Li in unit 10 



The chiral coefficients of relevance and electromagnetic pion mass difference are given 
in Tab. ^. Compared with Da Rold and Pomarol's case, the results do not significantly 
change much in our two cases. 



6. Conclusions 

In this paper, we considered holographic QCD beyond the leading order, by including two 
dim-6 dimension operators that go beyond the usual quadratic kinetic terms for the bulk 
gauge field Lm and Rm, and scalar field <I> 0]- We have studied the mass spectra. 
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decay constants of vector, axial and pseudoscalar sectors, and phenomenology of ai pir, 
p TTTT and oi TTj channels. In our work, we could achieve a non- vanishing branching 
ratio for ai ttj, which is a new feature compared with the usual holographic QCD in 
the leading order. We also calculated the electromagnetic form factor of a charged pion, 
(including the charge radius of a pion) which agrees with the experimental results up to 
~ 2 GeV'^. The numerical results are summarized in Table 1, and compared with the 
leading order results obtained by Da Rold and Pomarol Q denoted as the case A. We 
could achieve significant improvements in overall phenomenology of the tt — p — ai system 
by including the n and C terms. 

Let us remind ourselves that most studies based on the AdS/QCD approach are just 
the leading order calculations, starting from the bulk Lagrangian which is quadratic in the 
bulk gauge fields. Including the next-to-leading order corrections would be the next step to 
follow, and our present work makes such an attempt by considering dim-6 operators that 
reduce to the O(p^) operators after chiral symmetry breaking. Considering the improve- 
ment of overall phenomenology obtained in this work, it would be clearly desirable to have 
more systematic study of subleading corrections within AdS/QCD. 
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